Single particle spectra of charge transfer insulators by cluster perturbation theory - 

the correlated band structure of NiO 
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We propose a many-body method for band-structure calculations in strongly correlated electron 
systems and apply it to NiO. The method may be viewed as a translationally invariant version of 
the cluster method of Fujimori and Minami. Thereby the Coulomb interaction within the d-shells is 
treated by exact diagonalization and the d-shells then are coupled to a solid by an extension of the 
cluster perturbation theory (CPT) due to Senechal et al. The method is computationally no more 
demanding than a conventional band structure calculation and for NiO we find good agreement 
between the calculated single particle spectral function and the experimentally measured band 
structure. 

PACS numbers: 72.80.Ga,71.27.-a,79.60.-i 
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I. INTRODUCTION 

Band structure calculations based on the single par- 
ticle picture have enjoyed considerable success in solid 
state theory. Single particle picture here is meant to 
imply that the ground state is obtained by filling up ac- 
cording to the Pauli principle the energy levels calculated 
for a single electron in an 'effective potential'. The effec- 
tive potential thereby is usually constructed within the 
framework of the local density (LDA) or local spin den- 
sity approximation (LSDA) to density functional theory 
(DFT) 'I'l and despite the well-known fact that the eigen- 
values of the Kohn-Sham equations should not be identi- 
fied with the single-particle excitation energies of a sys- 
tem, the resulting band structures often give an almost 
quantitative description of angle-resolved photoemission 
spectroscopy (ARPES). 

However, there are also some classes of solids which defy 
such a description, most notably transition metal com- 
pounds with partially filled d and /-shells and strong 
Coulomb interaction between the electrons in these. A 
frequently cited example is NiO, where LSDA band struc- 
ture calculations correctly predict an antiferromagnetic 
and insulating ground state, but only a small 'Slater gap' 
of a fraction of an whereas experimentally NiO is 

an insulator with a bandgap of 4.3 eV|3| and stays so 
even above the magnetic ordering temperature. While 
DFT thus gives reasonable answers within its domain 
of validity - namely ground state properties - the non- 
correspondence between the Kohn-Sham eigenvalues and 
the single-particle excitation energies of the solid obvi- 
ously has to be taken literal for this compound (if the 
band gap is not read off from the LSDA band structure 
but expressed as the difference of ground state energies 
it is in fact possible to calculate it within the framework 
of DFT, as shown by Norman and Freeman 4] ). It is 
generally believed that the reason for the discrepancy is 
the strong Coulomb interaction between the electrons in 
the Ni Sd-shell, which leads to a substantial energy split- 
ting between configurations with different n. This 



leads to a very pronounced 'pinning' of the d-shell occu- 
pation number n, in the case of NiO to the value n = 8. 
Final states for photoemission or inverse photoemission 
then correspond to a single d-shell being in either a dJ , 
a or a configuration, in each case with very small 
admixture of configurations with other n. The corre- 
sponding 'defect' then may be thought of propagating 
through the crystal with definite k. This pinning of the 
electron number in both initial and final states cannot be 
reproduced by a wave function which takes the form of a 
simple Slater-determinant - such as the ground state de- 
duced from the Kohn-Sham equations. The situation is 
improved somewhat in the self-interaction corrected ver- 
sion of DFtIsIj'g'I, which renders a certain fraction of the 
d-orbitals completely localized, so that their occupation 
number in fact does become pinned - for the remain- 
ing delocalized d-orbitals, however, the problem remains. 
Another way to achive the pinning of the d-shell occu- 
pancy is the use of an orbital-dependent potential in the 
framework of the so called LDA-I-U method0,l8]. Speak- 
ing about gap values the calculations based on the GW- 
approximation'o', To*! also need to be mentioned - these 
give ah initio gap values which are in good agreement 
with experiment but do not seem to reproduce the high 
energy 'satellite' in the photoemisison spectrum. 
In addition to the pinning of the d-shell occupancy, the 
multiplet structure of the metal ion poses a problem for 
single-particle theories as well. It is quite well estab- 
lished that the multiplet structure (appropriately mod- 
ified by the crystal field splitting) of the isolated metal 
ion persists in the solids. Clear evidence for this point 
of view comes from the fact, that angle integrated va- 
lence band photoemission spectra^J ^ ^'^^ ^ X-ray 
absorption spectrap^ of many transition metal com- 
pounds can be reproduced in remarkable detail by config- 
uration interaction calculations solving exactly the prob- 
lem of a single d-shell hybridizing with a 'cage' of lig- 
ands. In these calculations it is crucial, however, that 
the intra-shell Coulomb repulsion is treated in full detail. 
While the cluster method is spectacularly successful for 



2 



angle-integrated quantities its 'impurity' character un- 
fortunately makes it impossible to extract the dispersion 
relations of fc-resolved single particle excitations. 
Actual dispersion relations in the presence of strong 
Coulomb interaction were first studied by Hubbard |13|, 
thereby taking an entirely different point of view as com- 
pared to the single particle picture on which conventional 
band structure calculations are based. Thereby the d- 
shells first are considered as isolated, and their affin- 
ity and ionization spectra obtained, thereby treating the 
Coulomb repulsion exactly. In his famous papers Hub- 
bard used a much simplified model, where the orbital 
degeneracy of the d-level was neglected whence ioniza- 
tion and afhnity spectrum of the 'half-filled' d-shell col- 
lapse to single peak each, with the two peaks separated 
by the Coulomb energy U . Upon coupling the individual 
atoms to the solid, the ionization and affinity states of 
the individual atoms then are systematically broadened 
to form the two 'Hubbard bands'. The coupling to the 
solid was achieved originally by the famous Hubbard 1 
approximation, but in fact this may be interpreted as a 
particularly simple form of the cluster perturbation the- 
ory (CPT), proposed by Senechal et a/.0], where the 
individual 'clusters' consist of just a single d-shell. This 
suggests immediately to relax Hubbard's simplifications 
and take into account the full complexity of a transition 
metal oxide including the orbital degeneracy of the d- 
shell, the full Coulomb interaction between c?-electrons 
in these and the sublattice of ligands. This is essentially 
the purpose of the present manuscript. 
An important complication is due to the sublattice of lig- 
ands. It has been shown by Fujimori and Minami[Tll| that 
in discarding altogether the sublattice of ligands (in the 
case of NiO: the oxygen atoms), Hubbard actually went 
one step to far in his simplification of the model. Namely 
Fujimori and Minami showed that the top of the valence 
band in NiO is composed of states, where a hole is pre- 
dominantly in an oxygen atom, but somehow 'associated' 
with an n-conserving excitation of a neighboring c?-shell - 
i.e. a magnon or a d — d exciton. This type of state might 
be viewed as generalization of a Zhang-Rice singlet 1 5j in 
the Cu02 planes of cuprate superconductors. It was then 
found by Zaanen, Sawatzky and Allen^^ that there is 
a crossover between this so-called charge-transfer insula- 
tor and a more conventional Mott-Hubbard insulator as 
a function of two key parameters, the Coulomb repulsion 
U between electrons in the d-shell and the charge transfer 
energy A, which are defined as 

£:(d" -> = U-A, 

£;(d"i^d"-i) = A. 

Strong experimental support for the picture proposed by 
Fujimori and Minami is provided by the resonance be- 
havior of the photoemission intensity as seen in photoe- 
mission with photon energies near the 2p — > Sd absoption 
threshold[l3lll[l3. 

Adopting this point of view and using a simplified 
'Kondo-Heisenberg' model, in which the charge degrees 



of freedom on Ni where projected out, Bala et al. |2Clf 
then obtained dispersion relations of quasiparticles in 
NiO which in fact do contain the key feature seen in 
ARPES|2lL |2^ : the coexistence of strongly dispersive 
oxygen bands on one hand and a complex of practically 
dispersionless (i.e.: massively renormalized) bands which 
form the top of the valence band on the other hand. 
In the present theory no reduction of the Hamiltonian 
to a t-J-type model is performed. Rather we use the 
same basic idea as in the treatment of the Kondo lat- 
tice in Ref. |2^: the system is divided into sub-units 
which are treated exactly and the hybridization between 
the sub-units is treated approximately. To do so, we de- 
fine the ground state for vanishing hybridization as the 
'vacuum state' and treat the charge fluctuations created 
by the hybridization as 'effective Fermions', for which 
an approximate Hamiltonian can be derived and solved. 
It has been shown in Ref. that the Hubbard-I ap- 
proximation for the single-band Hubbard model can be 
re-derived in this fashion if the sub-units are taken to be 
only a single site - including more complex 'composite 
particles' which extend over several unit cells then im- 
proves the agreement with numerical results. The gener- 
alized Zhang- Rice singlets discussed above may be viewed 
as such composite particles. It is shown in Appendix I of 
Ref. 25] that this treatment is in fact equivalent to the 
original cluster perturbation theory of Senechal et aL|l4j 
provided the sub-units into which the system is divided 
are non-overlapping. This last requirement poses a sub- 
stantial problem for transition metal oxides, because the 
rocksalt lattice of NiO cannot be easily divided into non- 
overlapping sub-units which are still amenable to exact 
diagonalization without artificially breaking a symmetry 
of the lattice (which would lead to artificial symmetry 
breaking in the band structure). We therefore need to 
adjust the concept of cluster perturbation theory to this 
situation - which also is an objective of the present work. 
The remainder of the paper is organized as follows: in 
section II we discuss a simplified ID model, in section HI 
we present the general theory, in section IV we apply the 
theory to the ID model and compare the obtained single 
particle spectra with result from exact diagonalization, 
in section IV we discuss the ARPES spectra of NiO and 
section V gives the conclusions. 



II. A SIMPLIFIED MODEL 

For a start we consider the following minimal version 
of a ID charge transfer model (see Figure ^ : 

H = -A^<,d,^, + C/^4^d,,^4,rf,^, 

+ II ( + ^-c- ) (1) 




FIG. 1: Schematic representation of the Hamiltonian 1^1 and 
its parameters. 

This model describes a ID chain consisting of strongly 
correlated metal((i) orbitals and uncorrelated ligand (p) 
orbitals. Henceforth we choose tpd as the unit of energy 
and unless otherwise stated set tpp — 0. The relevant 
filling (which we will consider henceforth) of the model is 
3 electrons (or 1 hole) per unit cell. While our goal ulti- 
mately is to study realistic models for compounds such as 
NiO, our motivation for studying this highly oversimpli- 
fied model is as follows: it is simple enough so that rea- 
sonably large clusters (up to 6 unit cells) can be treated 
by exact diagonalization (ED) and the obtained exact 
results for the single-particle spectral function then can 
serve as a benchmark for the analytical theory. The very 
simple nature of the model thereby is highly desirable, 
because it results in a small number of 'bands' so that 
the comparison with theory is more significant than e.g. 
in the case of NiO. 

The quantity of main interest is the photoemission and 
inverse photoemission spectrum, defined as 



<5(^ + (4"-i)-4")) 
5(c.-(4"+i)-i?("))) (2) 



where a € {xy, xz, yz, . . .} denotes the type of d-orbital, 
and ^[/"^ {E^^) denote the v^^ eigen state (eigen en- 
ergy) with n electrons - thereby v — Q corresponds to 
the ground state. The spectral function for p-electrons is 
defined in an analogous way. 

To get an idea how to construct an adequate theory it 
is useful to compare the paramagnetic mean-field solu- 
tion of the model (i.e. with A ^mf = A -|- {na)U) 
and the results of an ED calculation, see Figure |21 The 
mean-field solution gives two bands 
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The lower (fully occupied) one of these has predominant 
p-character, the upper (half-occupied) one has predom- 
inant d-character. In the spectra obtained by exact di- 
agonalization, the p-like bands persists with an almost 
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FIG. 2: Single particle spectral functions A''~\k,uj) and 
j4'+^(fc,a;) obtained by mean-field solution of the model and 
by exact diagonalization of a system with 6 unit cells. The 
wave vector k increases from the lowermost to the upper- 
most panel in steps of ^. 5- functions have been replaced by 
Lorentzians of width O.OStpd. The part to the left (right) of 
the vertical dashed line shows A'~'(fe,c<;) (A^'^\k,uj)). 



unchanged dispersion. This is not really surprising, be- 
cause an electron in the respective state moves predomi- 
nantly on the p-sublattice and thus will not feel the strong 
Coulomb interaction on the d-sites very much. The band 
with predominant d-character, on the other hand disap- 
pears completely in the exact spectra. There is a diffuse 
band somewhat below — A = —3tpd and a second one at 
the energy —A + U = 5tpd- Clearly, these two resemble 
the 'Hubbard bands' expected for a strongly correlated 
system and the respective final states have the character 
of an empty or doubly occupied d-orbital. In addition 
to these Hubbard bands, however, there is a third group 
of peaks at energies « +tpd which energetically is close 
to the p-like band, has a mixed p-d character and which 
does in fact form the first ionization states of the system. 
Its closeness to the p- band would seem to suggest that 
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the final states have a hole predominantly on the p-sites, 
but it also has a significant admixture of d-weight and 
moreover is closer to the Fermi energy than the p-like 
band. In the charge transfer system under consideration, 
the d-likc band thus actually splits up into three bands, 
rather than the two Hubbard bands which one might ex- 
pect. 

Despite the highly oversimplified nature of the ID 
model, there is actually already a clear analogy to NiO: 
LDA bandstructure calculations produce two well- 
separated band complexes, the lower one (i.e.: the one 
more distant from the Fermi energy) with predominant 
oxygen character, the upper one with Ni character. This 
is quite similar to the mean-field solution in Figure El 
The actual photoemission spectra, however, show first of 
all a broad structure at binding energies > 8eV below 
the top of the valence band. Resonant photoemission 
experiments 0, 0| show, that the final states observed 
in this energy range have predominantly d7 character - 
clearly they should be identified with the d-like band at 
—4:tpd in our model. Next, at binding energies —6eV and 
—AeV ARPES experiments^^! find a group of strongly 
dispersive bands which closely resemble the oxygen-like 
bands obtained from an LDA band structure calculation 

- they obviously correspond to the dispersive p-like 'rem- 
nant' of the free-electron band in our model. Finally, the 
top of the valence band in NiO is formed by a group of al- 
most dispersionless bands, whereby the mixed Ni (f and 
(PL_ character of these states is established by resonant 
photoemission JX 18] - these states then would corre- 
spond to the low intensity band which forms the top of 
the electron annihilation spectrum. The cluster calcu- 
lation of Fujimori and Minami 11] suggests, that these 
states should be viewed as hole-like 'compound objects' 
where a hole on oxygen is bound to an excited state of cP 

- a type of state that might be viewed as a generalization 
of the Zhang-Rice singlet in the Cu02 planes of cuprate 
superconductors. As already conjectured by Fujimori 
and Minami their 'compound nature' would make these 
quasiholes very heavy, which immediately would explain 
the lack of dispersion seen in the ARPES spectra 22] . 
Finally, inverse photoemission shows the presence of an 
upper Hubbard band in NiO which is also present in the 
spectra of the ID model. 

The above comparison shows, that in addition to p-like 
holes we will need three types of 'effective particles' to 
reproduce the correlated band structure of the model. 
The two standard Hubbard bands, which correspond to 
^n±i_jjjj.g gj-^al states are not sufficient here. To get an 
idea what these states should be, let us start from the 
ionic limit, tpp = tpd = 0. The ground state then cor- 
responds to a constant number of electrons, n, in each 
metal d-shell and completely filled ligand p-shells, see 
Figure Switching on the hybridization integral tpd 
then will produce charge fluctuations: in a first step, a 
hole is transferred into a p orbital, thus producing a 
state in d-orbital number z, see Figure|3lD. The state 
has an energy of C/ — A relative to the original c?" state. 
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FIG. 3: Charge fluctuation processes relative to the purely 
ionic configuration. 



and will become our first 'effective particle' - these 'par- 
ticles' form the unoccupied Hubbard band. In a second 
step, the p-like hole can be transferred into a c?-orbital 
i' ^ i, thus producing a state, see Figure Oh- The 
latter has an energy of +A relative to the d" state and 
provides the second type of 'effective particle' - actually 
the one that forms the 'satellite' in the spectral func- 
tion. Finally, the hole in i' can be transferred back into 
a neighboring p-level, thereby leaving the orbital i' in a 
state d"*, i.e. an eigenstate of d" other than the orginal 
one, see Figure Ell. The 'compound object' consisting of 
d"* and a hole in a neighboring p-orbital will be the third 
type of effective particle, its energy relative to the original 
d" state is « 0, i.e. appropriate to give top of the valence 
band. These states maight be viewed as generalizations 
of a Zhang- Rice singlet, or an extreme case of an (either 
spin- or orbital-like) 'Kondo object'. Level repulsion due 
to hybridization between the d"~^ and the d"*L states 
will push the latter up to higher energies relative to the 
p-like bands - in this way, these states become the first 
ionization states. In the following, we will try to give the 
above considerations a more solid theoretical foundation 
and apply them to the calculation of 'correlated band 
structures'. 



III. GENERAL THEORY 



We consider a typical transition metal oxide and re- 
strict our basis to the oxygen 2p-orbitals and transition 
metal d-orbitals. Taking the energy of the p-level as the 
zero of energy the single-particle terms in the Hamilto- 
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nian then take the form 

i.a,j.K, a 

Hid - (C' 4o,- + ^-c-), 

= E E (^".^ 4-,- + ^-^o, (3) 

where dj creates a spin-cr electron in the d-orbital 
a € {xy, xz, yz, . . . } on metal site i and ^ ^ creates 
an electron in p-orbital k G {x,y, 2} on oxygen site j. 
The Va,!3 combine charge-transfer energy and crystalline 
electric field. The Coulomb interaction between the d- 
electrons is 

E '^i:i:AAAA.^ (4) 

Cl,C2,C3,C4 

where we have suppressed the site label i-, C, — {m, a) and 
m G {—2. ..2} denotes the z-component of the orbital an- 
gular momentum. The matrix-elements VS^'S'^ can be 
expressed in terms of the 3 Racah-parameters, A B and 
C. 

For a start we take all the hybridization matrix elements 
t to be zero. In this limit each oxygen is in a 2p^ con- 
figuration, and each transition metal ion in one of the 
ground states of Hintra = + He with n electrons. In 
general, this ground state is degenerate, and we deal with 
this by choosing one of these ground states, which we call 
|$i^o)i for each metal ion - for example in the case of NiO 
we would choose the direction of the spin 5 = 1 of the 
d-shell to oscillate between the two sublattices, so as to 
describe the antiferromagnetic order in the system. We 
call |$i.o) the corresponding 'reference state' on transi- 
tion metal site i, it obeys Hintra\^i.o) = £'o"''|'&i,o)- In 
the following, we consider the product state of the |$i,o) 
and the completely filled oxygen-p sublattice as the 'vac- 
uum' of our theory. For the ID model |^ we introduce 
two d-sublattices A and B and choose 

1*^,0) = 4i' ^ 

so as to model the antiferromagnetic spin correlations in 
the system. 

Next, we assume that the hybridization between the sub- 
systems is switched on. This will create charge fluctua- 
tions in the vacuum: in a first step an electron from a 
p-shell will be transferred into one of the d-orbitals of the 
neighboring metal ion i, a process frequently denoted as 
— + Due to the many-body character of the 

Hamiltonian, the resulting state, dl^^\^i^o), in general 
is not an eigenstate of Hintra - rather, we can express it 



as a superposition of eigenstates: 

dUa\^.,o)=Y. C-^,o.,oAv). (5) 

v 

where |z/), v = 0,1, . . . v^ax are the eigenstates of the 
d-shell with n + 1 electrons. Since we are dealing with 
a single d-shell these can be obtained by exact diagonal- 
ization, which gives us the eigenenergies of the 

states Iv) as well as the coefficients Ci^ci,a,u- If the \v) 
are chosen to be eigenstates of Sz (as we will assume in 
all that follows) only a small fraction of the Ci^a,rT,v is 
different from zero. We now represent the state where 
the metal ion i is in the state \v) by the presence of a 
'book keeping Fermion', created by e|^^, at the site i. 
The spin index a thereby gives the difference in z-Spin 
between the state \v) and the reference state |$i,o) - in 
principle this is redundant, but we add it so as to make 
the analogy with a free-particle Hamiltonian more obvi- 
ous. An important technical point is, that in case the 
\v) are not eigenstates of 5*2 this labelling is not possi- 
ble - there may exist states \v) which can be reached by 
transferring an electron of either spin direction into the 
d-shell i. 

All in all, the charge fluctuation process then can be de- 
scribed by the Hamiltonian 

v;x. - E C (H4a..i*-o) 

a 

= E ^1" C!i,a,(7,iy- (6) 

a 

It is understood that states |z/) which have Ci^a,a,u — for 
both directions of a should be omitted from this Hamil- 
tonian. One can see, that the 'bare' hopping integrals tj'^ 
are multiplied by the coefficients Ci^a,a,u, which have a 
modulus < 1. The effective Fermions e|j^ thus in general 
have a weaker hybridization with the p-Orbitals than the 
original electrons, which will naturally lead to some kind 
of 'correlation narrowing' of all bands of appreciable d- 
character. 

In the ID model (Q) there is only one state with n + 1 
electrons, namely the state c?|^d]^jO) This has an energy 
of U — 2A, whence £i = U — A. For a site i on the A 
sublattice we thus identify el^\vac) = d|j^d-||0) and the 
part Hi becomes 

i7i = (f/-A)^ eli^a 
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plus an analogous term which describes the charge fluc- 
tuations on the sites of the B-sublattice. 
We proceed to the next type of state which is admixed 
by the hybridization. First, an electron from the metal- 
ion with ri + 1 electrons may be transferred back to the 
oxygen atom with the hole, i.e. d"^^L — > c?". If the 
metal ion i thereby returns to the reference state |$i,o) 
this process is described by the term 'iJ.c' in If the 
metal ion returns to an n-electron state |A) other than 
I'I'i.o) we should model this by a Bosonic excitation bj^. 
Here we will neglect these latter processes - this is pre- 
sumably the strongest and least justified approximation 
in the present theory. It implies for example that we are 
neglecting (in the case of NiO) the coupling to d—d exci- 
tons and the influence of the quantum spin fluctuations 
(spin waves) in the 3D S=l Heisenberg antiferromagnet 
formed by the Ni-moments. 

A second type of state can be generated by filling the 
hole in the oxygen shell j with an electron from a d-shell 
i' i, that means c?"L — > d"^^. This leaves the d-shell 
on i' in an eigenstate of n — 1 electrons, the net ef- 
fect is the transfer of an electron between the d-shells 
i' ~* z', i.e. precisely the process considered originally by 
Hubbard. We write 

d^,c.J'^>^,Q) =^ C^,a,aJ^^) (8) 

and model the shell i' being in the ionization state by 
the presence of a hole-like book-keeping Fermion, created 
by h\, ^ ^. In an analogous way as above we arrive at 
the following effective Hamiltonian to describe the second 
type of charge fiuctuation: 

a 

= E(0* '5.,.,.,.. (9) 

Since h\ ^ creates a hole-like particle, the presence of 
terms like h'^p^ is not unusual - these terms describe 
particle- hole correlations, not particle-particle correla- 
tions as in BCS theory. 

In the ID model, the only state with rt — 1 electrons is 
the empty site |0), which has the energy iJ = 0, whence 
e = A. We identify, for a site i on the B sublattice: 
h\^\vac) — |0) and the term H2 reads 

-^p'^E((?^Ui-?^lia)'^lT + ^-) (10) 



plus an analogous term describing the A sublattice. 
As already stated, the second type of charge transfer ex- 
citation describes the transfer of an electron between two 
d-shells via an intermediate state with a hole in an oxy- 
gen p-shell. If there are direct d-d transfer integrals, this 
process can also occur in one step. The respective part 
of the effective Hamiltonian can be constructed in an en- 
tirely analogous fashion as the parts above - since it is 
lengthy, we do not write it down in detail. 
We proceed to the last type of 'effective particle' that we 
will consider. If the d-shell on atom i is in an eigen state 
of n — 1 electrons (which would be described by the h\ ^- 
particle) it is possible that an electron from a neighboring 
p-shell is transferred to the d-shell, thereby leaving the 
d-shell i in an eigenstate |A) of n electrons other than the 
reference state |$i,o)- We will consider the 'compound 
object' consisting of the 'excited' d" state |A) on site i 
and a hole in a linear combination z/;^ ^ of p-orbitals on 
the nearest neighbors of i as a further effective particle, 
created by zj u- Here a S {xy, xz, yz, . . .} denotes the 
symmetry of the linear combination of p-orbitals, which 
is such as to hybridize with exactly one of the d] ^ ^ on 
site i. For an ideal tetrahedral cage of oxygen atoms 
around each metal ion there is exactly one combination 
^PI^ ^ for each a. The creation and annihilation of the 
z-particles can be described by the term 

i.fi.X.a 

f ^ - E''"^ - E''"^ - C 
tj,A,Q! — J^Q Set: 

Vi:,^o.,a = (MdlM (11) 

Here (a = (0|^q, ^[^^pp, "0,1, cr]|0) denotes the kinetic en- 
ergy of the combination ipj^ ^ - it can be expressed in 
terms of the integrals (ppcr) and (ppn). Also, = 
{0\di,a,a[Hpd,'4'l a a]\^) hybridization matrix ele- 

ment between i/ja a ^^id q a- ^'^d can be written in terms 
of (pdcr) and (pdn). 

In the ID model, the only 'excited' state on the A sub- 
lattice is the state rf|jO) The only possible combina- 
tion of p orbitals which hybridizes with a d-orbital is 

= 7f(( "^I-i<T )' ^^^"^^ " ~^PP- 

There is therefore just one z^^-like particle on the A- 
sublattice, which corresponds to site i being in the state 
J 0) and having an extra hole in the combination "01 1 o- ■ 
The total energy is i? = — A -I- tpp whence e — tpp and 
the corresponding Hamiltonian reads 

H3 = tpp J2 4^»T - ^(^/i,^ + H.c.) (12) 

and, again, a corresponding term for the B-sublattice. 
This concludes the types of state which we take into ac- 
count. We are thus assuming that the hole always is on 
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a nearest neighbor of the d"* state - this means that we 
truncate the 'Kondo cloud' which is not exactly true. In 
principle this approximation could be relaxed by includ- 
ing more complex composite particles but here we do not 
include these. 

In order for the mapping between the actual system and 
the 'book-keeping Fermions' to be a faithful one, we must 
require that the occupation of any d-site is either or 
1 - otherwise, the state of the respective d-shell is not 
unique. This implies that the book-keeping Fermions , 
and have to obey a hard-core constraint, in ex- 
actly the same way as e.g. the magnons in spin-wave 
theory for the Heisenberg antiferromagnet (HAF). It is 
well known, however, that linear spin-wave theory for 
the HAF, which neglects this hard-core constraint alto- 
gether and treats the magnons as free Bosons, gives an 
excellent description of the antiferromagnetic phase, even 
in the case of 5* = 1/2 and d = 2, where quantum fluc- 
tuations are strong. The reason is, that the density n of 
magnons/site obtained self-consistently from linear spin 
wave theory is still relatively small, whence the proba- 
bility that two magnons occupy the same site and thus 
violate the constraint is oc ^ 1. For the same reason 
we expect that relaxing the constraint in the present case 
and treating the book-keeping Fermions as free Fermions 
will be a very good approximation - its physical content 
is the assumption, that the probability of charge fluctu- 
ations is small, which is certainly justified in a Mott- or 
charge-transfer-insulator. For completeness we note that 
there is also a certain interference between the particle 
and the holes on oxygen in the sense that the respective 
creation and annihilation operators do not exactly anti- 
commute. Again, we neglect this, with the justification 
again being the very low density of the z^ and parti- 
cles. 

Adding the various terms in the Hamiltonian and the di- 
rect p — p hopping Hpp we obtain a Hamiltonian which 
describes the lowest order charge fluctuation processes 
while still being readily solvable by Fourier and Bogoli- 
ubov transform with the result: 

H=J2 ^k.c,a7,i,c,.7k,f,., (13) 

where ^ is a band index. Quantities of physical interest 
now can be readily calculated. Let us first discuss the 
electron count. We assume that the reference states for 
the d-shell have n electrons each. Then, the total number 
of electrons/unit cell is rie = n • + 6 • rio where nd(no) 
denote the number of metal (oxygen) atoms in the unit 
cell. On the other hand we have 

rie = n • rid + ^ ( l]pI,,-,«,aPfeJ,«,a 



Here the sums over j and i run over the p and li-shells 
in one unit cell and the equation follows readily from 
the electron/hole-like character of the various effective 
Fermions. This can be rewritten as 

1 V- t 

- y^(A<tot,(7 + hot,a)- (15) 
a 

where iJ,tot,iT denotes the total number of ionization states 
in the unit cell which can be reached by extracting a spin 
a electron from one of the |$i,o)) and and Xtot,iT the total 
number of 2:'^-particles in the unit cell which couple to 
one of these ionization states. If we assume that the band 
structure is spin independent (which is the case for NiO) 
these numbers must be independent of a and we obtain 
the following requirement for the electron number: 

E 7L,.7fe,a,a = N{6no + 2(Mtot + Xtot) ) • (16) 

fe,(T 

The number of occupied bands in the system thus is 
nocc = 3no -I- fitot + Atot. Since the total number of bands 
produced by our formalism is 'irio + vtot + IHot + ^tot , we 
find that the chemical potential falls exactly into the gap 
between the Srio + Utot + \ot bands which correspond 
(in the limit of vanishing hybridization) to the oxygen 2p 
states, the ionization states of the d-shell and the d^L- 
type states on one hand and the Vtot bands, which cor- 
respond to the affinity states of the rf-shells on the other 
hand. Obviously, this is the physically correct position 
of the chemical potential. 

The quantity of main interest to us, the single-particle 
spectral function Aa{k,Lxj) can be obtained from the 
eigenvectors of the Hamilton matrix once the resolution 
of the d-electron creation/annihilation operator is known. 
Were it not for the presence of the 2;'f-like particles, we 
could write 

(17) 

which is easily verified by taking matrix-elements be- 
tween the right an left-hand side. However, the presence 
of the z^-like 'particle' complicates this. Due to their 
'compound nature' the processes by which the electron 
annihilation operator couples to a ^'''-particle are rather 
complicated. For example, one might envisage a process 
in which an electron in a rf""*"^ configuration on site i is 
annihilated, leaving the d-shell in a d"* state (i.e. an 
eigenstate of other than the reference state |$i,o) on 
site i). Then, if simultaneously a hole happens to be 
present in a p-orbital next to site i, this process would 
create a z^-like particle on site i, leading to an opera- 
tor product of the type zjhitpia to describe this process. 
Similarly, if a d"* state is somehow created on a site i 
(this is not possible in the framework of the Hamiltonian 
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FIG. 4; Ground state energy and d-occupation as a function 
of U for A = 3, tpd = 1 as obtained by exact diagonalization 
of a system with 6 unit cells and from the present theory. The 
trivial contribution of —A has been subtracted off from the 
ground state energy. 

which we wrote down above - it would necessitate terms 
including the Bosonic excitations 6]^ discussed above) 
and a hole is created in a p-orbital next to this site, this 
would result in the creation of a z^-like particle on site 
i. This process could be described by a product of the 
type zj'i/iabi. Due to the fact that these processes all in- 
volve products of three operators one might expect that 
they lead predominantly to an incoherent continuum in 
the spectral function. All in all, we may thus expect 
that this type of process will not contribute substantially 
to the photoemission intensity for the dominant peaks. 
Clearly, the problem in calculating the spectral weight is 
a drawback of the theory - it should be noted, however, 
that there is a very clear physical reason for this problem, 
namely the 'compound nature' of the zj particles and this 
should be reflected in any theoretical description. 



IV. COMPARISON WITH EXACT 
DIAGONALIZATION 

Our theory involves a number of strong approximations 
which need to be checked in some way. Here we present a 
comparison of results obtained for the ID model and 
exact diagonalization of finite clusters. For the ID model 
systems with 6 unit cells easily can be solved exactly on 



a computer and we use these results as a benchmark to 
check our theory. The simplicity of the model actually 
makes the comparison more significant than e.g. in the 
case of a realistic model for NiO, because we expect only 
a small number of 'bands' whence any disagreement will 
be more obvious. 

Before we discuss results for physical quantities let us ad- 
dress one of the key approximations of our theory, namely 
the neglect of the hard-core constraint which in principle 
should be obeyed by the book-keeping Fermions. Solv- 
ing the ID model with A = 3, U — 6 gives the GS 
expectation values {d\d^) = 0.151, {h\h^) = 0.004 and 
{zlz^ = 0.0007. This implies that the probability for a 
violation of the constraint on any given d-site is « 0.03, 
i.e. entirely negligible. Simply relaxing the constraint 
thus is probably an excellent approximation. 
Next, Figure 01 compares the total energy/site and the 
d-occupancy, which is a measure for the charge-transfer 
form p ^ d, as obtained by exact diagonalization and 
from the theory. Obviously, there is good agreement. 
Figure |S1 shows a comparison between the single parti- 
cle spectral function obtained from the theory and by 
exact diagonalization of a system with 6 unit cells. To 
obtain a denser mesh of fc-points, spectra for a system 
with periodic and antiperiodic boundary conditions have 
been used for the exact diagonalization part, that means 
fc = 0, ^, ^ and TT have been calculated with PBC, the 
ones for A: = i-, -I and have been obtained with ABC. 

6 ' 2 6 

Although there is no rigorous proof for this, experience 
shows that combining spectra with PBC and ABC gives 
quite 'smooth' dispersion relations - as can also be seen 
in the present case. In order to suppress the Luttinger- 
liquid behaviour expected for ID systems, a staggered 
magnetic field of O.ltp^ was applied. The agreement be- 
tween theory and exact diagonalization then is obviously 
quite good. The dispersion and spectral character of the 
main 'bands' in the numerical spectra is reproduced quite 
well. The main difference concerns the very strong damp- 
ing of the lower Hubbard band at i? ~ —4:tpd, which ac- 
tually forms broad continuum rather than a well-defined 
band in the numerical spectra. Moreover, the upper Hub- 
bard band at « Atpd has some 'fine structure' in the 
numerical spectra, which is not reproduced by the the- 
ory. On the other hand, our theory does not include 
any damping mechanism such as the coupling to spin 
excitations, so one cannot expect it to reproduce such 
details. Another slight discrepancy concerns the band- 
width of the oxygen band at i? w 0, which is somewhat 
underestimated by theory. Apart from that and a few 
low-intensity peaks in the numerical spectra, however, 
there is a rather obvious one-to-one correspondence be- 
tween the bands in the theoretical spectra and the exact 
ones. Next, we consider the spectra for a nonvanishing 
p-p hopping, tpp = —1, see Figure El Again, there is 
good agreement between theory and numerics, with the 
main discrepancy being again the damping of the satellite 
and the fine structure of the upper Hubbard band. Still, 
there is a clear one-to-one correspondence between the- 
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FIG. 5: Single particle spectral functions A^~\k,uj) and 
j4''+'(fc,a;) obtained by the present theory and by exact di- 
agonalization of a system with 6 unit cells. The wave vector 
k increases from the lowermost to the uppermost panel in 
steps of -I, to that end the figure combines spectra obtained 
with periodic and antiperiodic boundary conditions. The part 
to the left (right) of the vertical dashed line shows A^~\k,Lij) 
{A(+\k,u;)). 



ory and exact spectra. An interesting check is provided 
by inverting the sign of tpp. One might expect at first 
sight that the only effect is to invert the dispersion of the 
p-like band. Inspection of the Hamiltonian H12() shows, 
however, that inverting the sign of tpp also affects the en- 
ergy of the z'I'-particle, and hence should lead to a shift 
of the corresponding band. The actual spectra in Figure 
13 then show, that this is indeed the case in the numerical 
spectra. The z-like band is shifted to higher energies by 
very nearly the amount of 2tpd expected from theory, so 
that the lowest hole-addition states now belong to the 
p-like band. The fact that the inversion of the sign of tpp 
has precisely the effect predicted by theory is a strong 
indication, that this is indeed the correct interpretation 
of the low energy peaks in the spectra. 
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V. THE BAND STRUCTURE OF NIO 

Summarizing the results of the preceeding section we 
may say that the theory reproduces the numerical spectra 
and the trends under a change of parameters remarkably 
well, an indication that despite its simplicity the the- 
ory really captures the essential physics of the two-band 
model. This is encouraging to apply it to a real mate- 
rial, NiO. In applying the above procedure to NiO we 
first performed a standard LDA band structure calcula- 
tion in the framework of the LMTO-method|23 for NiO 
(thereby assuming a paramagnetic ground state) and ob- 
tained the LCAO parameters by a fit. For simplicity no 
overlap integrals were taken into account. A comparison 
between the LDA band structure and the LCAO fit is 
shown in Figure |Hl(Jhe LDA result is essentially identi- 
cal to that of Ref.[22|), the hybridization integrals and 
site energies obtained by the fit are given in Table We 
have also obtained LCAO parameters for an antiferro- 
magnetic LSDA band structure, and those parameters 
which can be compared (such as the hybridization inte- 
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FIG. 8: LDA band structure for paramagnetic NiO and 
LCAO fit. The Fermi energy is taken to be zero. 



grals) do not differ significantly. All in all this procedure 
gives quite reliable estimates for the values of the various 
hopping integrals. The LDA bandstructure broadly can 
be divided into two complexes of bands: the lower one 
at energies between —8eV and —3eV has almost pure 
oxygen p-character. In other words, a hole in this bands 
would move almost exclusively on the oxygen sublattice 
and have only a very small probability to be on a Ni ion. 
One may thus expect, that these states persist essentially 
unchanged in the the correlated ground state. Next, the 
complex between —3eV and + leV has almost exclusively 
Ni 3d character. The LDA band structure thus would 
seem to suggest, that there are states, where a hole is 
moving essentially from one Ni site to another, which 
have a less negative binding energy - i.e. which are closer 
to the Fermi energy - than the states where the hole is 
moving in the oxygen sublattice. Clearly, in view of the 
value of the charge transfer energy A > 0, which is con- 
sistently suggested by a variety of methods 4,, Jl, 12], this 
is a quite wrong picture of the electronic structure. Next, 
we consider the Racah parameters B and C. These differ 
only slightly from their values for free ions and we took 
the values from Fujimori and Minami'll] B = 0.127 eV, 
C = 0.601 eV for and B = 0.138 eV, C = 0.676 eV 
for d^ . In general, these parameters are screened by co- 
valency between d-orbitals and ligandspTf but for sim- 
plicity we keep the 'bare' values. 

This leaves us with two parameters, which require a spe- 
cial treatment, namely the Racah parameter A, which is 
subject to substantial solid-state screening, and the dif- 
ference of site energies between the Ni 3d-level and the 
oxygen 2p-level. The Racah parameter A is related to the 
Coulomb energy U, which can be obtained from 'pure d- 



quantities' according to U ^ Eg+^ + - Here 

we used the values U — 8.7eV and A = 1.5eV. Sim- 
ilar values for U have been obatined by Fujimori and 
MinamifllJ from a cluster fit of the valence band pho- 
toemission spectrum, by van Elp et al. 12] from a clus- 
ter fit to the X-ray absorption spectrum and by Norman 
and Freemaui^] from density functional calculations. The 
value of A is somewhat small compared to others, which 
are around 2.5eV. 
Then, the problem of a single d-shell was solved by ex- 
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TABLE L Hybridization integrals and site-energies (in eV) 
obtained by a LCAO fit to the paramagnetic LDA band struc- 
ture of NiO. 
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FIG. 9: Single particle spectral densities (see l|5J for a def- 
inition) for antiferromagentic NiO obtained by the present 
theory. The momenta are along the (100) direction, the top 
of the valence band is the zero of energy. iS-functions are re- 
placed by Lorentzians of width O.OTSeV, the d-like spectral 
density is multiplied by a factor of 4. 



act diagonalization in the 7, 8 and 9-electron subspaces. 
The maximum dimension of the Hilbert space was 120 for 
n = 7. A nonvanishing CEF parameter lODq = Q.05eV 
was appUed in order to stabihze the correct t2geg ^^2g 
ground state for cfi in Oft,-symmetry. To account for the 
antiferromagnetic nature of the GS of NiO, we chose the 
reference state |$i,o) to be the Sz = I member of the 
^A2g multiplet on the Ni-sites of one sublattice, and the 
5"^ = — 1 member on the other one. Since we neglect spin- 
orbit couphng the direction of the spin quantization axis 
is arbitrary and has no influence on the spectral function. 
The kinetic energies of the t2g and Cg-like combinations 
of p-orbitals, which enter the energy of the z-like parti- 
cles are et^^ = (ppc) - (ppn) and = (ppn) - (ppa), the 
respective hybridization integrals are Tt^^ — 2{pdT:) and 
= V5{pda) . All in all, the rank of the effective Hamil- 
ton matrix to be diagonalized was « 250 i.e. quite mod- 
erate. To improve the agreement with experiment, the 
following minor adjustments of parameters were made: 
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FIG. 10; Comparison between the experimental peak disper- 
sions determined by ARPES in non-normal emission (taken 
from Figure 12 of Ref. p^D and the position of 'significant 
peaks' in the theroretical spectra. The labels on the 'bands' 
indicate a possible correspondence between experiment and 
theory. 



the p — p hybridization integrals were reduced by a factor 
of 0.8, and the d—p hybridization integrals were increased 
by a factor of 1.1. 

The full single particle spectral functions obtained along 
the (100) direction for antiferromagnetic NiO then is 
shown in Figure |^ It differs quite significantly from the 
what one would expect on the basis of the LDA band 
structure (see Figure JS))) but instead shows the same 
overall structure as in the ID model, compare Figure |21 
As was the case for the ID model one can broadly speak- 
ing distinguish four complexes of bands. At binding ener- 
gies < —8eV , there is a broad continuum of bands with 
strong d- weight. Analysis of the wave functions shows, 
that the respective states have (mainly) (i.e.: d^) char- 
acter, with some admixture of (i.e.: d^L) and (less) 
admixture of 02p character. Clearly these bands should 
be identified with the 'satellite' in the experimental NiO 
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spectra. By analogy with the ID model we may expect 
that these high energy states undergo substantial broad- 
ening as is indeed seen in experiment. In Figure El the 
satellite by and large disperses upwards as one goes away 
from r - Shen at aL[22| interpreted their data as show- 
ing a downward dispersion of the satellite. On the other 
hand this feature is rather broad and composed of many 
'subpeaks' so that it may be difficult to make really con- 
clusive statements about the dispersion of the spectral 
weight without a full calculation of the spectral weight, 
including the 'radiation characteristics' of the individual 
d-orbitals, final states effects etc. 

Next, there is a group of strongly dispersive bands of 
predominant 02p character, which closely resembles the 
lower complex of bands in the LDA calculation, see Fig- 
ure IHl In view of their almost pure oxygen character it 
is no surprise that these bands are hardly influenced by 
whatever happens on the Ni sites. Next comes a group 
of practically dispersionless bands which form the top of 
the valence band. They have mainly (i.e.: d^L) char- 
acter with some admixture of h'' (i.e.: cf). Due to their 
strong z^-character these bands probably are influenced 
most strongly by our approximation to omit any terms in- 
volving z^-operators in the spectral weight operator (|17|l . 
We may expect that taking the prcocesses discussed there 
will probably enhance the weight of these states and also 
add some more p-like weight to these peaks. 
The topmost peak is rather intense and actually com- 
posed of several 'subpeaks' - it is in fact the only feature 
in this energy range which shows signiflcant dispersion. 
Below this broad peak, there are several bands with lower 
intensity and practically no dispersion - all of this exactly 
as seen in the ARPES experiment by Shen et al.^^. Fig- 
ure 1101 shows a more detailed comparison of the disper- 
sion of 'significant peaks' in the photoemission part of 
the theoretical spectra with the experimental peak dis- 
persions as reported by Shen et ai. It can be seen, that 
the agreement is quite good. Along both (100) and (110) 
the main discrepancy is the position of bands C and D 
(or Di along (110)) which are somewhat higher in energy 
in the theory - still, the discrepancy is only a fraction of 
an eV. In view of the fact that we have used the simple- 
most set of parameters this is quite good agreement. The 
band portion Ea which is unusual due to its downward 
curvature has actually been observed by Shen at al. in 
normal emission (see Figure 6 of Ref. [22|). The part Eb 
seems to correspond to the experimental band E itself - 
it has rather low spectral weight for momenta close to F. 
Finally, FigurelTTI gives the dispersion of the 'sub-bands' 
of the broad structure A at the valence band edge. This 
fine structure has not yet been resolved experimentally 
as yet - however, Shen et al. found evidence for at least 
three 'subpeaks' and also for a quite substantial disper- 
sion, although this made itself felt only as a dispersion of 
the line shape of the broad peak. Looking at Figures 9 
and 10 of Ref. |22| it would appear that along (100) there 
is an overall 'upward' dispersion of the topmost peak A 
as one moves from T ~f X with two local maxima of the 



-0.2 - - -0.2 - 



-0.4 - , , ■ - -0.4 - 




Along (100} ^ ^ Along (110) 



FIG. 11: Fine structure of the broad peak which forms the 
top of the valence band structure in antiferromagnetic NiO. 
The symbols give the positions of peaks with appreciable 
weight. 



upper edge of A just after F und just before X with the 
whole band complex being most narrow approximately 
halfway between F and X. It can be seen already from 
Figure ^| that this dispersion of the peak-shape of A is 
reproduced quite well by theory. Similarly, along (110) 
the brod band complex seems to have its minimum width 
halfway between F and X. At least these qualitative re- 
sults are quite consistent with the dispersion in Figure 
im Clearly a more detailed study of the fine structure 
of feature A would provide an interesting check of the 
present and other theories for the band structure of NiO. 
Another stringent check for theory would be to unravel 
the orbital character of the individual fiat bands such as 
C and D by studying their intensity as a function of pho- 
ton polarization and energy. 

Finally, we mention the upper Hubbard band, with 
the corresponding final states having predominantly e^- 
character. The insulating gap has a magnitude of A.SeV, 
which is consistent with experiment Q. Figure IT?! shows 
the angle integrated (i.e.: k-integrated) photoemission 
and inverse photoemission spectrum. By and large there 
is reasonable agreement with experiment. The fact that 
theory puts the dispersionless bands C and D too close to 
the top of the valence bands leads to a too weak shoulder 
on the negative binding energy side of the 'main peak' at 
the top of the valence bands. 



VI. CONCLUSION 

In summary, we have presented a theory for the single 
particle-excitations of charge transfer insulators. The ba- 
sic idea is to interpret the charge fluctuations out of the 
purely ionic configuration as 'effective Fermions' and de- 
rive and solve an effective Hamiltonian for these. This is 
the same physical idea which is underlying both the Hub- 
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FIG. 12; Momentum integrated spectral weight for antifer- 
romagnetic NiO. To simulate a photoemisison spectrum the 
Lorentzian broadening has been taken energy dependent ac- 
cording to 5 = OAeV + {uj ~ leV) * 0.1. 



bard I approximation and the cluster perturbation theory 
and, as demonstrated above, when applied to a realis- 
tic model of a charge-transfer insulator these methods, 
which so far have been restricted to more 'model-type' 
systems, give indeed quite satisfactory agreement with 
experiment. The key approximation, namely to treat 
the Hubbard-like operators describing the charge fluctu- 
ations as free Fermions thereby is well justified, because 
of the low density of these effective Fermions, which ren- 
ders their (strong) interaction largely irrelevant. A sys- 
tematical way to relax this approximation would be the 
T-matrix approach, as demonstrated by Kotov et al\2^. 
It should be noted that the calculation is computationally 
no more demanding than a conventional band structure 
calculation and can be 'automated' almost completely. 
The weakest link in the chain thereby is the necessity to 
perform an LCAO-fit to an LDA band structure. 
One important conceptual problem is the necessity to 
break the symmetry which originates from the degen- 
eracy of the ground state multiplet of a single transi- 
tion metal ion and choose the 'reference states' |$i,o) 'by 
hand'. However, one might as well consider choosing an 
ansatz for these reference states which takes the form of 
a linear combination |<i>i,o) = OLiii\v), where the sum 
extends over the GS multiplet, and determine the coeffi- 
cients tti^ from the requirement of minimum total energy. 
In this way spin and orbital ordering could be studied in 
much the same way as lattice parameters are optimized 
in conventional LDA calculations and since all 'ingredi- 
ents' for the Goodenough-Kanamori rules'2^ are taken 
into account, this may be a quite promising method. 
Since spin-orbit coupling also can be trivially included 
in the exact diagonalization of the isolated d-shells one 
might even hope to address magnetic anisotropics and/or 
anisotropic exchange interactions. A procedure for the 
improvement of CPT calculations on model Hamiltoni- 



ans which is similar in spirit has been proposed by Pot- 
thoff et a/. 113]. 

One major drawback of the theory clearly is the approx- 
imate nature of the calculation of the spectral weight. 
It should be noted, however, that there is a very clear 
physical reason for this problem, namely the 'compound 
nature' of the ZRS-like states which form the top of the 
valence band. If the present interpretation of these states 
is the correct one, basically any theory will face similar 
problems. One possible way out would be to derive a 
version of the original CPT which can work with site- 
sharing clusters. 

Finally, we would like to discuss the relationship between 
our theory and previous workers in the field. Manghi et 
aL[3]| and Takahashi and Igarashi 32] have calculated 
the quasiparticle band structure of NiO more along the 
lines of conventional many-body theory. Starting from a 
paramagnetic LDA band structure (Ref. i31|]) or an an- 
tiferromagnetic Hartree-Fock band structure (Ref. 32]) 
these authors added a self-energy constructed within the 
local approximation to three-body scattering theory. The 
obtained band structures show the same 'large scale fea- 
tures' as the one obtained here, but there are also signif- 
icant differences, particularly so near the top of the va- 
lence band. More detailed comparison with experiment 
seems necessary to discuss the merits of the various the- 
ories. 

Next, there is a clear analogy between the present theory 
and the cluster method of Fujimori and Minami^^ and 
va Elp at a/. ^3]- With the exception of cPl7 states in 
the photoemissions spectrum the present theory employs 
the same type of basis states as the cluster calculations. 
The only difference is, that we designate one of the de- 
generate ground states of d" as a 'vacuum state' and 
interpret the other states as 'deviations' from this vac- 
uum state. Those deviations which carry the quantum 
number of an electron then are considered as effective 
free Fermions. As discussed above, the low density of 
these effective Fermions probably make this a very good 
approximation. 

There is also an obvious relationship between the present 
theory and the work of Unger and Fuldefs^- Us- 
ing t he p rojection technique developed by Becker and 
FuldeyJ] these authors constructed an equation of mo- 
tion for single-particle spectral functions of the Cu02 
plane, which is very similar to the ones which would be 
obtained from our effective Hamiltonians. 
Finally we address the work of Bala et al. 20], which is 
very similar in spirit to the present theory. These au- 
thors derived a 'Kondo-Heisenberg'-like model operating 
in the subspace of (PL_ type states by eliminating - via 
canconical transformation - the charge fluctuations be- 
tween states of the type A2g<P)L and states of the 
type d'' (their theory was concerned with the motion of 
a single hole in an 02p orbital). Accordingly, their the- 
ory produced (in addition to the free-electron-like 02p 
bands) two weakly dispersive bands - one for each of 
the 'flavours' Cg and t^g whereby the flavour stands for 
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the symmetry of the linear combination of 02p orbitals 
around a given Ni site. Thereby Bala at al. actually 
went one step beyond the present theory by taking into 
account the coupling of 02p-like holes to the antiferro- 
magnetic magnons - which is omitted in the present the- 
ory. In the present theory, no canonical transformation is 
performed, so that also the high energy features (satellite 
and upper Hubbard band) are reproduced. Moreover, we 
also take the excited multiplets of and their covalent 



mixing with the <f multiplets into account, whence we 
obtain a larger number of ZRS-like bands - consistent 
with experiment. Experimentally the impact of the cou- 
pling to magnons which is ignored in the present theory 
but treated accurately in the work of Bala at a. could 
be studied only by considering the 'fine structure' of the 
broad peak at the valence band top. These states seem 
to have an appreciable dispersion which might or might 
not be influenced by the coupling to magnons. 
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